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Abstract 

The far-infrared absorption of a two-dimensional electron gas with a square- 
lattice modulation in a perpendicular constant magnetic field is calculated 
self-consistently within the Hartree approximation. For strong modulation 
and short period we obtain intra- and intersubband magnetoplasmon modes 
reflecting the subbands of the Hofstadter butterfly in two or more Landau 
bands. The character of the absorption and the correlation of the peaks to 
the number of flux quanta through each unit cell of the periodic potential 
depends strongly on the location of the chemical potential with respect to the 
subbands, or what is the same, on the density of electrons in the system. 
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Indications of the Hofstadter subband structure [||,|| of the Landau bands of a two- 
dimensional electron gas (2DEG) in square modulated lateral superlattices have been be- 
lieved to be found directly ^ and indirectly 0] in transport measurements. This opens the 
important question of whether the subband structure can directly be detected in far-infrared 
(FIR) absorption measurements. Calculations of the ground state properties of the laterally 
modulated 2DEG indicate that the strong screening effects would hamper an observation of 
even the coarse structure of the Hofstadter butterfly in, e.g. magnetocapacitance or magne- 
toresistance measurements, unless very short periods and strong modulation would be used 
1^. But what about FIR measurements? 

The FIR-absorption is determined by the self-consistent field which describes the dynam- 
ical screening of the incident external field. Since, as calculations |^ show, screening is also 
very important in the ground state, screening effects should be treated on equal footing in 
both the ground state and in the FIR response. That an inconsistent treatment of screening 
effects may lead to qualitatively wrong results in a 2DEG with strong lateral modulation 
is known from the case of quantum dots. There, according to the generalized Kohn's the- 
orem [§,0, Coulomb interaction effects on the ground state and dynamical response of a 
parabolically confined electron system cancel, so that the dot responds at the bare single- 
electron frequencies. This experimentally confirmed fact has been reproduced only by those 
many-body calculations which have treated the interaction effects in the ground state and 
in the FIR response on the same footing. Results of such calculations have been published 
for isolated quantum dots and for unidirectionally modulated 2DEG |]TO| , pT| , but to our 
knowledge no corresponding results are available for the square modulation.^ 

The square lateral superlattice with period L is spanned by the lattice vectors R = 
mil + ^12? where li = Lie, and I2 = Ly are the primitive translations of the Bravais lattice 
B; n,m & Z. The reciprocal lattice 71 is spanned by G = Gigi + G'2g2, with gi = 27ry//i, 
g2 = 27rx/Z2, and Gi,G2 G Z. The ground-state properties of the interacting 2DEG in 
a perpendicular homogeneous magnetic field B = Bz and the periodic potential V{r) = 
Vo{cos((7ix) + cos{g2y)} are calculated in the Hartree approximation. The Hartree single 
electron states \a) and their energies Ea are labelled by the quantum numbers {n^,/!, z/} = 
a, where ri; > is a Landau band index, /i = (6'i + 27mi)/p, v = {62 + 27rn2)/g, with 
rii G /i = {0, • ■ - p — 1}, 722 G /2 = {0, ■ ■ ■ g — 1}) G [— 7r,7r], and pq E N is the number of 
magnetic flux quanta through the lattice unit cell. The Hartree states are expanded in terms 
of the symmetric basis functions constructed by Ferrari and used by Silberbauer 
and the present authors 0. The Ferrari basis is complete provided that (/x, z/) 7^ (vr, vr) for 
all (ni,n2) G /i x I2. The Hartree potential "felt" by each electron and caused by the total 
electronic charge density of the 2DEG, —ens{r), and the positive neutralizing background 
charge, +enb, together with the external periodic potential V{r) do not couple states at 
different points in the quasi-Brillouin zone 9 = {61,62) G {[— vr, tt] x [— 7r,7r]}. However, 
states at a given point depend, in a self-consistent way, on states in the whole zone through 
ns(r) ||T3|j5|]. The periodic potential broadens the Landau levels into Landau bands that 



The collective excitations of a 2DEG with a two-dimensional magnetic- field modulation have 
been studied by X. Wu and S. E. Ulloa in Phys. Rev. B 47, 10028 (1993), but without considering 
interaction effects in the ground state. 



2 



due to the comensurability conditions between L and the magnetic length / = (ch/eBY^'^ , 
are spht into pq subbands. The resulting energy spectrum, the Hofstadter butterfly [|I],0, 
retains essentially it's complicated gap structure but is strongly reduced in symmetry by the 
electron-electron interaction |^ and coupling to higher Landau bands caused by the periodic 
potential |T^. The symmetry depends on the mean electron density, or, equivalently, the 
chemical potential. 

In order to calculate the FIR absorption of the 2DEG we perturb this by a monochromatic 
external electric fleld: 

Eext(r, t) = ~^^0 ||^^ Q| *^^P {^(^ + G) • r - iujt}. (1) 

Here we do not restrict the dispersion relation for the external fleld, C(j(k + G), to that of 
a free propagating electromagnetic wave, instead we allow for the more general situation in 
which the external fleld is produced as in a near-fleld spectroscopy or in a Raman scattering 
set-up. The power absorption is found from the Joule heating of the self-consistent electric 



fleld with = 4>ext + 4>ind, 

P(k + G, a;) = [|k + G|0.e(k + G, uj)<j>Uk + G, u)] . (2) 

The induced potential 4>ind is caused by the density variation 6ns{r) due to (psc, which can 
then in turn be related to the external fleld by the dielectric matrix 

^eG,G'(k,cu)0,c(k+G',cu) = 0exi(k + G, cj). (3) 

G' 

The dielectric tensor, eG,G'(k, t^^) = 5g,g' — K|k+G| ^G,G'(^^i is determined by the suscep- 
tibility of the 2DEG 

XG,G'(k,w) = ■^^Id9j:a,a' fli'^^i^^) J^^'^^'^^ (M + G) , . 

where k is in the flrst Brillouin zone, k is the dielectric constant of the surrounding medium, 

in which /° is the equilibrium Fermi distribution, i] ^ 0^, and 

^'J'(k) = («'(^')|e-^'-1«W)- (6) 

Special care must be taken with respect to the symmetry of the the wave functions corre- 
sponding to the Hartree states \a{0)) when translating them across the boundaries of the 
quasi-Brillouin zones. The peaks in the absorption spectra represent absorption of collective 
modes of the 2DEG. Single-electron transitions and collective oscillations can be identifled 
as zeros of the real part of the determinant of the dielectric matrix. For the collective os- 
cillations the imaginary part simultaneously vanishes but becomes very large for the single 
electron transitions indicating the strong damping of such processes. Here we prefer to 
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calculate P(k + G,uj) for a monochromatic external field instead of searching for zeros of 
det{eG,G'} in order to avoid the complication of multiple branches of the magnetoplasmon 
dispersion for different reciprocal lattice vectors G. The response of the system can also be 
analyzed in more detail by imposing a certain symmetry or polarization on the external field. 
Although we study collective oscillations here, it is still possible to identify in many cases the 
most important single-electron transitions contributing to the collective mode, either from 
the dispersion or by restricting the type of single electron transitions in the calculation. 

In the numerical calulations we use GaAs parameters, m* = 0.067mo, and k = 12.4. The 
power dissipation of the 2DEG is made possible by retaining a small but finite imaginary 
part for the frequency in the susceptibility, Eqs. (^) and (|^). We choose r] = uJbO. In 
order to examine the fine structure of the magnetoplasmon dispersion we have chosen a 
monochromatic external perturbing field (|1]) with G = and k in the first Brillouin zone. 

The absorption of a homogeneous 2DEG (Vq = 0) is compared in Fig. [l| with the first- 
order dispersion in k of a magnetoplasmon, u"^ = u"^ + 2Tie'^nsk/ {nm*). Kohn's theorem 
manifests itself by the fact that only one peak is visible for k [0, and for finite wave 
vectors k the magnetoplasmon interacts with harmonics of the cylcotron resonance, resulting 



in so the called Bernstein modes 16-18 



To make the Hofstadter subband structure of the Landau bands discernible, we have to 
resort to short periods L = 50 nm and strong modulation Vq = 4 meV, just as we had to in 
order to observe the subbands in the thermodynamic density of states ^j. In addition, the 
structures are very sensitive to the location of the chemical potential with respect to the 
subbands, i.e. to the filling factor v of the Landau bands. Figure |]a shows the absorption 
for u = 1/2, and the low frequency part is reproduced in Fig. In Fig. ^ the subband 
dispersion in one quasi- Brillouin zone along several values of 6*2 is projected onto the 6'i-axis, 
thus indicating the bandwidths and the position of the chemical potential. Since one and 
a half subband in the lowest Landau band are filled, and it is not inhibited by a selection 
rule, we can see three intraband magnetoplasmons. The fc-dispersion allows the following 
identification in Fig. ^d: The lowest frequency peak is the plasmon with contributions from 
transitions between the upper two subbands, the center peak gets contributions from the 
two lowest subbands and the high frequency peak has contributions connecting the lowest 
and the top subband. Since only two subbands of the lowest Landau band are occupied (of 
which one is partially filled), only two peaks are found for the interband magnetoplasmon 
with frequency just above Uc- The absorption for the same situation but for a different filling 
factor, i.e. u = 5/6, is seen in Fig. ^ Now all the three subbands of the lowest Landau band 
are at least partially occupied and three interband magnetoplasmon absorption peaks are 
clearly discernible. For u = 1 the low frequency intraband peaks vanish and the interband 
peak with the lowest frequency accounts for almost all of the oscillator strength. The 
structure of the Hofstadter spectrum can be seen for other values of pq in the calculations. 
In Fig. ^ the absorption for the cases pq = 2 and = 4 is seen for quite different parameters. 
The double structure expected for pq = 2 is visible for relatively large period, L = 100 nm, 
and weak modulation Vq = 0.4 meV, while two of the four peaks for the case of pq = 4 are 
just appearing for Vq = 3 meV, and L = 50 nm. This reflects the stronger screening of the 
2DEG observed in the ground state calculations ^ for higher values of pq, where a strong 
modulation and short period were needed to see structure in the thermodynamic density 
of states caused by the Hofstadter subbands. In addition, here the filling factor u not only 
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determines the strength of the screening and thus the relevance of selection rules, but it is 
also essential to control in order to occupy subbands such that single electron transitions 
can take place contributing to the collective oscillations. 

It is interesting to note that in the case of a nonvanishing potential modulation, Vq 7^ 0, 
the absorption P{w) has a nontrivial structure, even for the case k — * 0. The condition 
for the validity of the generlized Kohn's theorem can thus be violated by either spatially 
modulating the time-dependent external electric field (0) or by the static potential ^(r) the 
2DEG resides in. 

In summary, the FIR-absorption of a 2DEG in a short-period strongly modulated square 
lateral superlattice shows the underlying Hofstadter subband structure. The structure can 
be seen in the absorption due to either intra- or interband magnetoplasmons by tuning the 
density or the filling factor within the lowest Landau band. The intraband magnetoplasmon 
peaks in the absorption have weak oscillator strengths compared to the interband absorption 
and are situated in a region on the energy scale that is difficult to detect with present day 
FIR technology. 
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FIGURES 



FIG. 1. The power absorption P{k\,uj) for a homogeneous 2DEG (Vq = 0). The dashed straight 
Unes represent the frequencies of the cyclotron resonances u) = nujc and the dashed curve is the 
dispersion of the magnetoplasmon in first order with respect to k\. L = 200 nm, pq = 1, u = 1, 
m* = 0.067mo, T = 1 K, k = 12.4, and hujc = 0.178 meV. 

FIG. 2. (a) The power absorption P(/ci,u;) (in arbitrary units) for a modulated 2DEG 
(Vb = 4 meV). The solid curve in the (/ci, a;)-plane represents the dispersion of the magneto- 
plasmon of a homogeneous 2DEG in first order with respect to ki. (b) The low frequency region 
repeated from (a). The left peak decreases with increasing wave vector, ki. (c) The dispersion of 
the two lowest Landau bands in the quasi-Brillouin zone along several values of 02 projected on 
the 6i axis (each band is split into three subbands). The chemical potential is indicated by a solid 
horizontal line, pq = 3, L = 50 nm, u = 1/2, and Huc = 8.57 meV. 

FIG. 3. (a) The power absoprtion P(/ci,u;) for a modulated 2DEG {Vq = 4 meV). The solid 
curve in the {ki,Lo)-p\ane represents the dispersion of the magnetoplasmon in a homogeneous 2DEG 
in the first order with respect to ki. (b) The dispersion of the two lowest Landau bands in the 
quasi-Brillouin zone along several values of 02 projected on the 0i axis (each band is split into 
three subbands). The chemical potential is indicated by a solid horizontal line, pq = 3, L = 50 nm, 
u = 5/6, and huc = 8.57 meV. 

FIG. 4. The power absorption P{ki,u:) for a modulated 2DEG in two cases; (sohd) pq = 2, 
L = 100 nm, u = 1/2, huc = 1.43 meV, and Vq = 0.4 meV, and (dashed) pq = 4, L = 50 nm, 
u = 7/8, huc = 11.4 meV, and Vq = 3 meV. kiL = 0.6, and T = 1 K in both cases. 
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